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ABSTRACT 

In this second paper of a series, we present a treatment procedure for data cubes 
obtained with the Spectrograph for Integral Field Observations in the Near Infrared 
of the Very Large Telescope. We verified that the treatment procedure improves signif¬ 
icantly the quality of the images of the data cubes, allowing a more detailed analysis. 
The images of the Bry and H 2 A 212 I 8 emission lines from the treated data cube of the 
nuclear region of NGC 5643 reveal the existence of ionized and molecular-gas clouds 
around the nucleus, which cannot be seen clearly in the images from the non-treated 
data cube of this galaxy. The ionized-gas clouds represent the narrow-line region, in 
the form of a bicone. We observe a good correspondence between the positions of the 
ionized-gas clouds in the Bry image and in an [O III] image, obtained with the Hubble 
Space Telescope , of the nuclear region of this galaxy convolved with an estimate of the 
point-spread function of the data cube of NGC 5643. The morphologies of the ionized 
and molecular gas seem to be compatible with the existence of a molecular torus/disc 
that collimates the active galactic nucleus (AGN) emission. The molecular gas may 
also flow along this torus/disc, feeding the AGN. This scenario is compatible with the 
unified model for AGNs. 

Key words: techniques: imaging spectroscopy - galaxies: individual: NGC 5643 - 
galaxies: Seyfert - galaxies: nuclei 


1 INTRODUCTION 

Data cubes are data sets with two spatial dimensions and 
one spectral dimension. Data cubes enable the analysis of 
spectra at different spatial positions of a given object, or 
to perform cuts along the spectral direction, and thus to 
produce images of the object at different wavelengths. The 
analysis of data cubes in astronomy in the last decades has 
increased significantly with the construction of more instru¬ 
ments like Integral Field Spectrographs and Fabry-Perot 
spectrographs. Although the wealth of information in this 
kind of data is very significant, the performed analysis is usu¬ 
ally affected by the existence of artefacts, like high spatial- 
frequency noise or ‘instrumental fingerprints’ (see Fig. 5 of 
Neumayer et al. 2007). Therefore, a proper removal of these 
artefacts is recommended if one wishes to perform a cleaner 
analysis. 

The Spectrograph for Integral Field Observations in the 
* E-mail: robertobm@astro.iag.usp.br 


Near Infrared (SINFONI; Eisenhauer et al. 2003; Bonnet et 
al. 2004) of the Very Large Telescope ( VLT ) provides 3D 
imaging spectroscopy in a spectral range of 1.05 — 2.45pm 
(covering the spectral bands J, H and K). This instru¬ 
ment, actually, consists of a combination of the Adaptive 
Optics Module, developed by the European Southern Obser¬ 
vatory, with the Spectrometer for Infrared Faint Field Imag¬ 
ing (SPIFFI), developed by the Max-Planck-Institute for 
extraterrestrial Physics. SINFONI has three different fore¬ 
optics, which result in field of views (FOVs) of 8.0 arcsec 
x 8.0 arcsec (hereafter 8.0 arcsec), 3.2 arcsec x 3.2 arc- 
sec (hereafter 3.2 arcsec) and 0.8 arcsec x 0.8 arcsec (here¬ 
after 0.8 arcsec). These three fore-optics have spatial pixels 
(spaxels) of 0.125 arcsec x 0.250 arcsec (spaxel scale of 250 
mas), 0.05 arcsec x 0.10 arcsec (spaxel scale of 100 mas) and 
0.0125 arcsec x 0.0250 arcsec (spaxel scale of 25 mas), re¬ 
spectively. The fore-optics with an FOV of 0.8 arcsec allows 
observations with a spatial resolution of about 0.056 arcsec 
in a wavelength of 2.2pm (which is a near-diffraction limited 
resolution), when the adaptive optics (AO) is applied. The 
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Table 1. Data used throughout the paper to present the treatment procedure. 


Target 

Programme ID 

Name of PI 

Observation date 

Spaxel scale (mas) 

Filter 

Exposure time (s) 

Number 
of frames 

HIP 044245 

60.A-9235(A) 

- 

2009 April 8 

250 

K 

2 

i 

HIP 048480 

074,B-9012(A) 

Eisenhauer 

2005 March 16 

100 

K 

30 

i 

HIP 049220 

074.B-9012(A) 

Eisenhauer 

2004 December 22 

25 

K 

10 

i 

HIP 096155 

60.A-9235(A) 

- 

2007 August 14 

250 

H+K 

1.5 

i 

HIP 091481 

60.A-9235(A) 

- 

2006 September 17 

100 

H+K 

30 

i 

HIP 024761 

60.A-9235(A) 

- 

2011 August 6 

25 

H+K 

2 

i 

NGC 5128 

074.A-9011(A) 

Eisenhauer 

2005 March 24 

100 

K 

900 

9 

NGC 5128 

075.B-0490(A) 

Cappellari 

2005 April 2 

250 

J 

300 

4 

NGC 5128 

075.B-0490(A) 

Cappellari 

2005 April 2 

250 

H 

300 

4 

NGC 5128 

075.B-0490(A) 

Cappellari 

2005 April 2 

250 

K 

300 

4 

NGC 5643 

083.B-0332(A) 

Hicks 

2009 April 8 

250 

K 

600 

5 


SINFONI detector is composed by one CCD with 2048 x 
2048 pixels. 

The incident light in the FOV of SINFONI passes 
through the pre-optics region, which has three main sub- 
areas. In the first of these sub-areas, the background ther¬ 
mal emission is suppressed. The second sub-area is a mo¬ 
torized filter wheel that selects one of the four broad-band 
filters available: J, H, K and H+K. Since the H+K filter 
corresponds to the H and K bands together, it provides the 
largest spectral coverage, but the lowest spectral resolution. 
The third sub-area is the motorized optics wheel, which se¬ 
lects one of three lenses corresponding to the available fore¬ 
optics (with FOVs of 8.0, 3.2 and 0.8 arcsec). After passing 
through the pre-optics region, the light is transferred to the 
image slicer, where, at first, it is divided into several parts 
(the slitlets) by a group of 32 small mirrors (the small slicer). 
After this division, the light passes by a second group of big 
mirrors (the big slicer), which reconfigure all the slitlets in 
a way that, at the end of this process, they are all disposed 
in a shape similar to a slit. The light is then transferred 
to the collimator, where it is redirected to the region con¬ 
taining the diffraction gratings (which generate the spectra). 
Finally, the light is transferred to the detector of the instru¬ 
ment. Each of the 32 slices is projected on to 64 detector 
pixels, explaining why the spaxels on the sky are rectangu¬ 
lar. 

This is our second paper focused on describing data 
treatment procedures for data cubes obtained with differ¬ 
ent instruments. In Menezes et al. (2014, hereafter Paper 
I), we presented a procedure for treating data cubes (after 
the data reduction) obtained with the Near-Infrared Inte¬ 
gral Field Spectrograph (NIFS; McGregor et al. 2002). In 
this paper, we present a similar procedure for treating SIN¬ 
FONI data cubes. Although the procedures described here 
and in Paper I have essentially the same steps, there are 
many differences in the way each one of these steps should 
be applied to NIFS and to SINFONI data cubes, in order 
to obtain optimized results. All the steps in this sequence 
were applied using scripts written in Interactive Data Lan¬ 
guage! (idl). The data used in the paper to illustrate the 
treatment procedure are described in Section 2. In Section 
3, we present details about the correction of the differential 

1 http://www.exelisvis.com/ProductsServices/IDL.aspx 


atmospheric refraction (DAR). In Section 4, we describe the 
spatial re-sampling of the data. The Butterworth spatial fil¬ 
tering, to remove high frequency components from the im¬ 
ages of the data cubes, is discussed in Section 5. The process 
of ‘instrumental fingerprint’ removal is presented in Section 
6. The Richardson-Lucy deconvolution, to improve the spa¬ 
tial resolution of the data cubes, is described in Section 7. 
The analysis of the data cubes of NGC 5643 (used as a sci¬ 
entific example to show the benefits provided by the data 
treatment) is discussed in Section 8. Finally, in Section 9, 
we present a summary of the procedure and our conclusions. 


2 OBSERVATIONS AND DATA REDUCTION 

To illustrate most of the treatment procedure, we used ob¬ 
servations of the standard stars HIP 044245, HIP 048480 and 
HIP 049220, available in the SINFONI public access data 
archive^- All the observations were taken in the K band, 
with different fore-optics. 

For the analysis of the spatial displacements of the ob¬ 
jects along the spectral axis of SINFONI data cubes (Section 
3), we used observations of the standard stars HIP 096155, 
HIP 091481 and HIP 024761, which were taken in the H+K 
band, with different fore-optics. We chose observations of 
stars in this specific spectral band because it has the highest 
spectral coverage, which makes it easier to analyse the spa¬ 
tial displacements along the spectral axis of the data cubes. 
To show the process for ‘instrumental fingerprint’ removal 
(Section 6), we used observations of the nuclear region of the 
galaxy NGC 5128 (Centaurus A), in the J, H and K bands 
in different fore-optics. 

Finally, as a scientific example illustrating the benefits 
obtained with the data treatment (Section 8), we analysed 
observations of the nuclear region of the galaxy NGC 5643, 
in the K band. The AO was applied to all the observations 
taken with the fore-optics with FOVs of 3.2 and 0.8 arcsec; 
however, it was not applied to the observations taken with 
the fore-optics with FOV of 8.0 arcsec. The details of all the 
data used throughout this paper are shown in Table [T] 

Standard calibration data of linearity lamp, distortion 


2 http://archive.eso.org/wdb/wdb/eso/sinfoni/form 
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Figure 1. Average spectra and images of intermediate-wavelength intervals of the data cubes of (top) HIP 049220, (middle) HIP 048480 
and (bottom) HIP 044245, obtained after data reduction, using the standard pipeline. Note that the contours of the stars are not well 
defined in the images. In addition, high spatial-frequency components can also be observed. 


fibre (used in the data reduction to compute spatial dis¬ 
tortions and to perform a spatial rectification), flat lamp, 
arc lamp and sky field were obtained from the public data 
archive. Calibration spectra of standard stars were also re¬ 
trieved for the reduction of the data cubes of NGC 5128 
and NGC 5643. The data reduction was performed with 
the GASGANO softwai'43 and included the following steps: 
correction of bad pixels, flat-field correction, spatial recti¬ 
fication (including correction for spatial distortions), wave¬ 
length calibration, sky subtraction and data cube construc- 


3 http://www.eso.org/sci/software/gasgano.html 


tion. One data cube was obtained for each one of the follow¬ 
ing standard stars: HIP 044245, HIP 048480, HIP 049220, 
HIP 096155, HIP 091481 and HIP 024761. The data cubes 
of HIP 044245 and HIP 096155 were reduced with spaxels 
of 0.125 arcsec x 0.125 arcsec (hereafter 125 mas), the data 
cubes of HIP 048480 and HIP 091481 were reduced with 
spaxels of 0.05 arcsec x 0.05 arcsec (hereafter 50 mas) and 
the data cubes of HIP 049220 and HIP 024761 were reduced 
with spaxels of 0.0125 arcsec x 0.0125 arcsec (hereafter 12.5 
mas). Throughout this paper, we will refer mainly to the 
original spaxel scales (25, 100 and 250 mas) of the data 
cubes of these standard stars (not to their HIP identifiers), 
as this is the most important factor in our discussion. Nine 
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Figure 2. Graphs of the coordinates X c and Y c of the centres of the stars (top) HIP 024761, (middle) HIP 091481 and (bottom) HIP 
096155 as a function of the wavelength. The values of A x and A y represent the spatial shifts along the horizontal and vertical axis, 
respectively, from the positions of the stars at the b lue extremities of the data cubes. The red curves are third degree polynomials, the 
green curves were obtained with the equations from lBdnsch fe Potulskl d 1998tl and the blue curves were obtained with the task ‘ refro ’ 
from the ‘SLALIB’ package. The green and blue curves are almost coincident; therefore, it is very difficult to differentiate them. 


data cubes were obtained for NGC 5128, with spaxels of 50 
mas, and five data cubes for NGC 5643, with spaxels of 125 
mas. Finally, all data cubes were flux calibrated and telluric 
corrected within iraiQ, using appropriate standard stars (as¬ 
suming they have blackbody spectra). Fig. [T] shows images 
of intermediate-wavelength intervals of the data cubes of 
HIP 049220, HIP 048480 and HIP 044245, obtained after 
data reduction, using the standard pipeline, as well as the 
corresponding average spectra. 


4 http://iraf.noao.edu 


3 CORRECTION OF THE DAR 

The first step in our data treatment sequence is the correc¬ 
tion of the DAR effect, which causes a displacement of the 
observed object along the spectral axis of a data cube. A 
detailed description of the DAR is given in Paper I. More 
details can also b e found in iBonsch fc Potulskil d 1998tl and 
iFilippenkol (1 19821 ). 

Before this procedure is discussed here, however, it is 
necessary to evaluate the significance of the DAR effect in 
SINFONI data cubes. This is particularly important in this 
case because, since the DAR is stronger in the optical than 
in the infrared, this effect is usually ignored in data cubes 
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Figure 3. RGB composite images of the data cubes of the 
stars (top) HIP 024761, (middle) HIP 091481 and (bottom) HIP 
096155. In all images, the blue colour corresponds to the wave¬ 
length range 1.4615 - 1.4685 /im, the green colour corresponds 
to the wavelength range 1.7920 - 1.7990 /im and the red colour 
corresponds to the wavelength range 2.4485 - 2.4555 /im. 


obtained in the infrared. In order to evaluate how significant 
is the spatial displacement of structures along the spectral 
axis of SINFONI data cubes, we determined the coordinates 
X c and Y c of the centres of the stars HIP 024761, HIP 091481 
and HIP 096155, which were observed in the H+K band, 
with spaxel scales of 25, 100 and 250 mas, respectively. Fig[5] 
shows the graphs of the coordinates X c and Y c of the centres 
of these stars as a function of wavelength. 

An analysis of Fig. [2] reveals that the total spatial dis¬ 
placements, along the H+K band, of the data cubes with 
spaxel scales of 25, 100 and 250 mas are ~ 100, ~190 and 


~270 mas, respectively, along the horizontal axis, and ~ 110, 
~40 and ~50 mas, respectively, along the vertical axis. The 
total displacements of the star observed with spaxel scale 
of 25 mas along the horizontal and vertical axis are almost 
twice the predicted spatial resolution (56 mas, with the use 
of the AO) for the corresponding fore-optics (with an FOV 
of 0.8 arcsec). This indicates that these displacements are 
relevant for studies involving the H and the K bands to¬ 
gether. However, the spatial displacements are comparable 
to the predicted spatial resolution even if only one of these 
spectral bands is used; therefore, we conclude that the DAR. 
effect may be significant for SINFONI data cubes obtained 
with the fore-optics with an FOV of 0.8 arcsec. In the case 
of the star observed with spaxel scale of 100 mas, the spa¬ 
tial displacement along the horizontal axis was also almost 
twice the typical spatial resolution (100 mas, with the use of 
the AO) obtained with the corresponding fore-optics (with 
an FOV of 3.2 arcsec); therefore, such displacement is rele¬ 
vant for studies involving the H and the K bands together. 
However, similarly to what was observed with the star with 
spaxel scale of 25 mas, even the spatial displacement along 
one of these two spectral bands is still comparable to the 
typical spatial resolution of this fore-optics. Based on these 
observations, we conclude that the DAR effect may be also 
significant for SINFONI data cubes obtained with the fore¬ 
optics with an FOV of 3.2 arcsec. Finally, in the case of the 
star observed with spaxel scale of 250 mas, it is important 
to remember that the AO was not applied during the ob¬ 
servation (this is usually the case when the corresponding 
fore-optics is used). Therefore, the predicted spatial resolu¬ 
tion for this fore-optics depends fundamentally on the seeing 
conditions of the observations. Considering that, we can see 
that the observed spatial displacements of this data cube 
along the horizontal and vertical axis are smaller than the 
typical spatial resolutions obtained with this fore-optics. As 
a consequence, we conclude that the correction of the DAR 
effect may not be necessary in SINFONI data cubes obtained 
with the fore-optics with an FOV of 8.0 arcsec. 

The entire discussion above indicates that the correc¬ 
tion of the DAR may be relevant in some cases, particularly 
when the fore-optics with FOVs of 0.8 and 3.2 arcsec are 
used. In Fig. [3l we can see RGB composite images showing 
different wavelength intervals of the data cubes with spaxel 
scales of 25, 100 and 250 mas. The scientific necessity of the 
correction of the DAR in the data cubes with scales of 25 
and 100 mas is evident in these RGBs. On the other hand, 
we can also see that the correction of the DAR may not be 
really necessary in the data cube with scale of 250 mas. 

Fig. [5] revealed that the DAR may be important even 
when only one spectral band is analysed. However, the spa¬ 
tial displacements caused by this effect are obviously more 
significant when more than one band is taken into account. 
A study that often involves the entire continuum of the J, H 
and K bands together is the spectral synthesis. Analysis of 
spectral lines may also be compromised by the DAR if the 
purpose is to compare images of lines located in different 
spectral regions. This is the case in the analysis of emission 
line ratios, for example. 

The data cubes with spaxel scales of 25, 100 and 250 
mas were observed with considerably high zenith distances: 
62°08, 63°28 and 63°56, respectively, corresponding to air- 
masses of 2.14, 2.22 and 2.25, respectively. The DAR will 
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Figure 4. Graphs of the coordinates Y c of the centres of the 
stars (top) HIP 024761, (middle) HIP 091481 and (bottom) HIP 
096155, keeping the original size of the spaxels along the verti¬ 
cal axis (25 mas for HIP 024761, 100 mas for HIP 091481 and 
250 mas for HIP 096155). The values of A y represent the spatial 
shifts along the vertical axis from the positions of the stars at 
the blue extremities of the data cubes. The red curves are third 
degree polyno mials and the green c urves were obtained with the 
equations from IBonsch fc Potul ski] (1998). 


be less significant for observations with lower zenith dis¬ 
tances and, in such cases, even studies using observations 
taken with the fore-optics with FOVs of 0.8 and 3.2 arcsec 
or studies involving the continua in different spectral bands 
simultaneously may not require any correction. 

Similarly to what we discussed in Paper I for NIFS data 
cubes, the significant impact of the DAR in some SINFONI 


data cubes may be somewhat unexpected, as this effect is 
usually ignored in infrared 3D spectroscopy. However, one 
should keep in mind that AO provides very high spatial res¬ 
olutions to SINFONI data cubes. Therefore, although the 
DAR is indeed smaller in the infrared than in the optical, 
this effect may be relevant in SINFONI data cubes due to 
the high spatial resolutions. The DAR will have more im¬ 
pact on future data cubes obtained with the next generation 
of extremely large telescopes (of 25-39 m apertures), as the 
AO-corrected point-spread functions (PSFs) in these obser¬ 
vations will be considerably smaller. 

The script we use for the DAR correction shifts each 
image of the data cube in order to remove any spatial dis¬ 
placement of the structures along the spectral axis. This cor¬ 
rection can be achieved using two strategies: the theoretical 
approach, in which the spatial di splacements are calculated 
usi ng the theoretical e quations of lBonsch fe Potulskil (Il998l) 
and lFilippenkol l|l982l ). and the practical approach, in which 
the spatial displacements are measured directly. These two 

strategies are described in further detail in Paper I. _ 

_ The theoretical equations from IBonsch fc Potulskil 

(1 19981) assume a plane-parallel atmosphere. However, a more 
precise calculation is performed by the task ‘refro’ from the 
‘SLALIB’ package, which takes into account the curvature 
of the Earth’s atmosphere and an approximation to the 
change in atmospheric temperature and pressure with al¬ 
titude. Fig. [2] shows third degree polynomials fitted to the 
point s, the theoretical c u rves o btained with the equations 
from IBonsch fe Potulskil (Il998l ) and the theoretical curves 
obtained with the task ‘refro’ from the ‘slalib’ package (im¬ 
plemented in Fortran environment). 

Fig. [2] reveals no significant difference between the 
theoretical curves obta ined with the equations from 
IBonsch fe Potulskil (Il998l ) and with the task 1 refro ’ from the 
‘slalib’ package. These theoretical curves also do not repro¬ 
duce exactly the values of X c and Y c of the centres of the 
stars along the spectral axis. The highest discrepancies be¬ 
tween the observed values and the theoretical curves for the 
data cubes with spaxel scales of 25, 100 and 250 mas are ~5, 
~55 and ~130 mas, respectively, along the horizontal axis, 
and ~11, ~16 and ~20 mas, respectively, along the verti¬ 
cal axis. These discrepancies are very common in SINFONI 
data cubes. Similar results were also obtained in Paper I for 
NIFS data cubes. Based on that, we conclude that the the¬ 
oretical curves (considering a plane-parallel atmosphere or 
even taking into account the curvature of the atmosphere) 
are not the most precise way to reproduce the spatial dis¬ 
placements of structures along the spectral axis of SINFONI 
data cubes. As a consequence, the practical approach has a 
higher precision than the theoretical approach for remov¬ 
ing the DAR from SINFONI data cubes. Nevertheless, one 
should keep in mind that the observed discrepancies men¬ 
tioned before are all smaller than the predicted spatial reso¬ 
lutions of the observations in the corresponding fore-optics. 
We have never seen discrepancies much higher than these. 
Therefore, we can say that the theoretical approach is also 
a sufficiently precise way for removing the DAR from SIN¬ 
FONI data cubes. In some cases, there may be real physical 
effects that could cause the spatial displacement of an ob¬ 
ject along the spectral axis of a data cube. In such cases, the 
practical approach should not be used. 

It is important to try to determine the cause of the dis- 
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Figure 5. Graphs of the coordinates X c and Y c of the centres of the stars (top) HIP 024761, (middle) HIP 091481 and (bottom) HIP 
096155 as a function of the wavelength, after the correction of the DAR, using the practical approach. The values of A,„ and A y represent 
the spatial shifts along the horizontal and vertical axis, respectively, from the expected positions of the stars, after the correction of the 
DAR. 


crepancies between the observed values of X c and Y c and the 
theoretical curves of the DAR effect, detected in Fig. [2] They 
were not caused by an imprecise calculation of the DAR, as 
even the detailed calculation performed by the task ‘ refro’ 
did not reproduce the observed spatial displacements. As ex¬ 
plained in Paper I, a possible cause for the discrepancies is 
the spatial re-sampling performed during the data reduction. 
Indeed, the data cubes obtained with the fore-optics with 
FOVs of 0.8, 3.2 and 8.0 arcsec were reduced with spaxels 
of 0.0125 arcsec x 0.0125 arcsec (hereafter 12.5 mas), 0.05 
arcsec x 0.05 arcsec (hereafter 50 mas) and 0.125 arcsec x 


0.125 arcsec (hereafter 125 mas), respectively; however, the 
original size of their spaxels were 0.0125 arcsec x 0.025 arc¬ 
sec, 0.05 arcsec x 0.1 arcsec and 0.125 arcsec x 0.25 arcsec, 
respectively. Therefore, the data reduction did not alter the 
size of the spaxels along the horizontal axis but reduced by 
half their sizes along the vertical axis. To check if this spa¬ 
tial re-sampling affected significantly the values of Y c of the 
centres of the standard stars analysed here, we performed a 
new reduction of these data cubes, but keeping the size of the 
spaxels along the vertical axis unchanged. Then, we deter¬ 
mined again the values of Y c of the centres of the stars along 
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the spectral axis of the data cubes. Fig. [4] shows the results 
and it is easy to see that the previously observed discrepan¬ 
cies between the observed values of Y c and the theoretical 
curves remained essentially the same. 

Based on Figs. 2 and 4, we can say that the discrep¬ 
ancies between the observed values of X c and Y c and the 
theoretical curves are not due to an imprecise calculation of 
the DAR effect or to the spatial re-sampling applied during 
the data reduction. We were unable to determine the exact 
cause of these discrepancies. We believe that an instrumen¬ 
tal cause for this behaviour is a possibility. The fact that 
the observed discrepancies appear to increase from the fore¬ 
optics with an FOV of 0.8 arcsec to the fore-optics with an 
FOV of 8.0 arcsec for stars observed with very similar zenith 
distances seems to support the idea of an instrumental effect. 
A proof of this hypothesis, however, is beyond the scope of 
this work, due to the necessity of the analysis of many data 
cubes, obtained with different observational parameters. 

We applied the correction of the DAR, using the prac¬ 
tical approach, to the data cubes of the standard stars with 
spaxel scales of 25 mas, 100 mas and 250 mas. The graphs 
with the values of X c and Y c of the centres of the stars in 
the obtained data cubes, as a function of the wavelength, 
are shown in Fig.[o] It is easy to see that the procedure was 
very effective in removing the DAR. Since the expected val¬ 
ues, after the correction, for X c and Y c , respectively, were 

31.5 and 24.5 for the data cube with spaxel scale of 25 mas, 

43.5 and 16.5 for the data cube with spaxel scale of 100 mas 
and 32.0 and 30.0 for the data cube with spaxel scale of 250 
mas, we can say that the precision (1 <t) of our algorithm 
was higher than 3 mas in all the corrections performed here. 
Our tests showed that precisions between 1 and 10 mas are 
expected from our algorithm. 

We also applied our algorithm to the data cubes of 
HIP 049220, HIP 048480 and HIP 044245, which were ob¬ 
served in the K band, with spaxel scales of 25, 100 and 250 
mas, respectively. These data cubes are used to describe the 
methodologies in the next sections. As mentioned in Section 
2, only one exposure for each one of these stars was retrieved. 
However, when there is more than one exposure per object, 
the best results are obtained when these exposures are com¬ 
bined in the form of a median, after the correction of the 
DAR. 


4 SPATIAL RE-SAMPLING 

The second step in our sequence for the treatment of SIN- 
FONI data cubes is the spatial re-sampling. We use an al¬ 
gorithm that reduces the size of the spaxels in the images of 
the data cubes, conserving the surface flux. Our algorithm 
is based on the task ‘frebin’ of the idl Astronomy User’s 
Librarj0. However, other tasks may also be used, as long as 
the surface flux is conserved. After the spatial re-sampling, 
we perform an interpolation of the values in each one of the 
lines and columns of the re-sampled images. Besides provid¬ 
ing a better visualization of the contours of the structures in 
the images, the main improvement achieved with this proce¬ 
dure is related to the final spatial resolution, which is higher 

5 http://idlastro.gsfc.nasa.gov 


when the Richardson-Lucy deconvolution (the last step in 
our sequence - Section 7) is applied to re-sampled images. 
It should be noted that the spatial re-sampling of an image 
does not change its spatial resolution, but only improves its 
appearance. 

In Paper I, we give details about the necessity of the 
use of an interpolation after the re-sampling of an image and 
about the fact that better results are obtained by applying 
the spatial re-sampling procedure, instead of simply reduc¬ 
ing the data cubes with smaller spaxels. Since this previous 
discussion is entirely analogous to SINFONI data cubes, it 
will not be repeated here. 

We usually apply the spatial re-sampling (followed by 
an interpolation) to SINFONI data cubes (observed with 
any fore-optics) in order to double the number of spaxels 
along the horizontal and vertical axis. Therefore, if data 
cubes observed with the fore-optics with scales of 25, 100 
and 250 mas are reduced with spaxels with sizes of 12.5, 50 
and 125 mas, respectively, then, the final sizes of their spax¬ 
els, after the spatial re-sampling, are 6.25, 25 and 62.5 mas, 
respectively. Fig.[S]shows images of intermediate-wavelength 
intervals of the data cubes with spaxel scales of 25, 100 and 
250 mas, before and after the spatial re-sampling, following 
the strategy mentioned above. 

In Fig, [6] we can see that the spatial re-sampling im¬ 
proved the appearance of the images of the data cubes. How¬ 
ever, we can also see the existence of some high spatial- 
frequency components, which appear as horizontal stripes, 
in the images of the re-sampled data cubes. A part of these 
high-frequency structures may have an instrumental origin 
and, therefore, was possibly already present in the data 
cubes before the spatial re-sampling. However, it is prob¬ 
able that some of them were introduced by the spatial re¬ 
sampling procedure. It is known that a spatial re-sampling 
introduces high spatial-frequency components in the images 
and that could explain, at least partially, the horizontal 
stripes detected in the images of the re-sampled data cubes. 
This side effect of the spatial re-sampling process is not 
problematic, as the high spatial-frequency structures are ef¬ 
ficiently eliminated by the Butterworth spatial filtering (Sec¬ 
tion 5). 

To analyse in further detail the high-frequency compo¬ 
nents introduced by the spatial re-sampling, we calculated 
the Fourier transforms of all the images in Fig. [6] The mod¬ 
uli of the obtained Fourier transforms are shown in Fig. [7] 
The bright elliptical regions farther from the centre of the 
Fourier transforms represent the high spatial-frequency com¬ 
ponents introduced by the spatial re-sampling. The defini¬ 
tion of the Nyquist frequency we are using is based on the 
spaxel sampling of each image. One interesting aspect of 
Fig- El is that the high-frequency components of the images 
are much more intense along the vertical axis than along the 
horizontal axis. This characteristic can also be seen in the 
re-sampled images in Fig. [S] as the high-frequency compo¬ 
nents appear mainly as horizontal stripes and not vertical 
stripes in the images. This behaviour is very common in 
SINFONI data cubes and represents a difference from what 
we have observed in NIFS data cubes (Paper I), in which the 
high-frequency components along the horizontal and verti¬ 
cal axis have similar intensities. Another important feature 
of the Fourier transforms in Fig. [3 is that the region cor¬ 
responding to the low spatial frequencies has an elliptical 
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Figure 6. Images of intermediate-wavelength intervals of the data cubes of the stars (top) HIP 049220, (middle) HIP 048480 and 
(bottom) HIP 044245, (left) before and (right) after the spatial re-sampling. 


shape. Again, this is characteristic of SINFONI data cubes 
and represents a difference from NIFS data cubes, in which 
the low-frequency region in the Fourier transforms has a 
shape similar to the product of an ellipse by a rectangle. 

Another possible side effect of the spatial re-sampling of 
the images in a data cube is the degradation of the signal- 
to-noise (S/N) rati o of the spec tra. However, in our pre¬ 
vious experiences (M enezesl I 2 OI 2 I) . we have never observed 
measurable considerable degradations. Therefore, consider¬ 
ing the benefits provided by the spatial re-sampling and the 
non-significant side effects, we conclude that the spatial re¬ 
sampling of SINFONI data cubes is recommended. 


5 BUTTERWORTH SPATIAL FILTERING 

The third step in our sequence for the treatment of 
SINFONI data cubes is the Butterworth spatial filtering 
(lOonzalez fc Woodj|2003 ). which is performed in the fre¬ 
quency domain. This procedure is described in detail in Pa¬ 
per I. We use a script that applies the filtering to all the 
images of the data cube, in order to remove the high spatial- 
frequency noise. 

The general formula of a low-pass Butterworth filter 
of order n is given by equation (7) of Paper I. Since the 
low-frequency region of the Fourier transforms of SINFONI 
data cubes have typical elliptical shapes (see Fig. 0, the 
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Figure 7. Moduli of the Fourier transforms of the images in Fig. [fi] Note that the high-frequency components introduced by the spatial 
re-sampling appear as bright regions farther from the centre of the Fourier transforms of the re-sampled data cubes. 


best choice for the Butterworth filter is an elliptical one, 
which is given by 


1 




1 + 

/( u-u n \ 2 i l «-»n \ 2 


V v a ) 1 V b ) 



In the previous equation, the point (uo,vo) is located 
at the centre of the Fourier transform, n is the order of the 
filtering, a is the cut-off frequency along the horizontal axis 
and b is the cut-off frequency along the vertical axis. Our 


tests revealed that, considering the eccentricities of the But¬ 
terworth filters involved, it is very difficult to detect mea¬ 
surable differences between the results obtained using an 
elliptical filter and using a simple circular filter (which is 
equivalent to assume a = b in equation 1). Therefore, we 
use a circular filter throughout this paper. 

Fig. [7] indicates that a cut-off frequency of approxi¬ 
mately 0.21 Ny seems to be adequate for the Butterworth 
spatial filtering of the data cube of the standard star ob¬ 
served with spaxel scale of 100 mas. Following the same 
approach we u sed in Paper I, we applie d a spatial wavelet 
decomposition dStarck fc Murtaghll200d l. using the A Trous 
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Figure 8. Images of the wavelet data cubes (a) Wo, (b) Wi, (c) W 2 , (d) W 3 , (e) W 4 and (f) Wc of HIP 048480 (with spaxel scale of 
100 mas), in an intermediate-wavelength interval. 


algorithm, to each image of this data cube. For more details 
about the wavelet decomposition of data cubes and about 
the A Trous algorithm, see Paper I. Fig. 0 shows images of 
an intermediate-wavelength interval of the obtained wavelet 
data cubes Wo, Wi, W 2 , W 3 , W 4 and Wc- The high spatial- 
frequency components appear very clearly in the image of 
Wo as horizontal stripes. The interference of these struc¬ 
tures in the image is so significant that the star cannot be 
seen directly. Part of these high-frequency components was 
probably introduced by the spatial re-sampling, although it 
is also possible that at least some of these stripes have an 
instrumental origin (as discussed in Section 4). In fact, the 
horizontal stripes were also easily detected in the image of an 
intermediate-wavelength interval of this data cube, after the 
spatial re-sampling (see Fig. O- The high spatial-frequency 
structures are much less abundant in the image of Wi and 
are practically absent in the other images in Fig. [H] 

We then applied the Butterworth spatial filtering to the 
Wo data cube of the standard star observed with spaxel scale 
of 100 mas, using a filter given by equation (1) with a cut¬ 
off frequency of 0.21 Ny and with n = 2, n = 3, n = 4, 
n — 5 and n — 6. Fig. [9] shows images of an intermediate- 
wavelength interval of Wo, filtered with these values of n. 

A simple analysis of Fig. [9] indicates that all the fil¬ 
tering procedures applied to Wo, with different values of 
n, removed a substantial amount of high spatial-frequency 
structures and, as a consequence, the star became visible in 


all filtered images of Wo- However, the images of the data 
cubes filtered with orders higher than 2 show the existence 
of ‘concentric rings’ around the central star. It is also easy 
to see that the higher the value of n, the more intense and 
more numerous are the rings. As discussed in Paper I, this 
behaviour can be explained by the fact that, based on the 
convolution theorem, a Butterworth spatial filtering of an 
image is equivalent to the convolution of the image with the 
inverse Fourier transform of the Butterworth filter. Consid¬ 
ering this and the aspect of the inverse Fourier transforms 
of the Butterworth filters with different values of n (shown 
in Fig. EB- the reason for the existence of concentric rings 
in the images of data cubes filtered with orders higher than 
2 becomes very clear. 

There are some points that are worth mentioning here. 
First, although the concentric rings introduced by the But¬ 
terworth spatial filtering with high orders were evident in 
Fig. [9] it is much more difficult to detect them in images of 
original SINFONI data cubes (without a wavelet decompo¬ 
sition), because they are usually immersed in the PSFs of 
the observations. They were so clearly seen here because we 
used the Wq data cube in this discussion. Besides that, con¬ 
sidering the way these rings are introduced in the images, 
they are only detected around point-like sources. Finally, in 
a Butterworth spatial filtering, there is an issue associated 
with the periodicity of Fourier transforms. Such periodicity 
is a consequence of the way a Fourier transform is defined. 
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Figure 9. Images of an intermediate-wavelength interval of (a) the non-filtered Wo data cube of HIP 048480 (with spaxel scale of 100 
mas) and of the corresponding filtered data cubes, using (b) n = 2, (c) n = 3, (d) n = 4, (e) n = 5 and (f) n = 6. 


The results obtained by a convolution between two func¬ 
tions with and without periodicity are different. Therefore, 
since the Butterworth spatial filtering of an image is equiv¬ 
alent to the convolution between the image and the inverse 
Fourier transform of the filter, we conclude that this process 
is affected by the periodicity of the Fourier transforms. To 
avoid this problem, we apply a ‘padding’ procedure, which 
will not be discussed here, but a de tailed description can be 
found in lGonzalez fe Woodsl J2002I) . 

Since the Butterworth spatial filtering with n = 2 re¬ 
moved a substantial amount of high spatial-frequency noise 
and did not introduce visible artefacts in the images, we can 
say that this seems to be the most adequate value of n to be 
used in the Butterworth spatial filtering of SINFONI data 
cubes. However, a careful analysis of Fig. [9] reveals that, al¬ 
though the Butterworth spatial filtering of Wo with n = 2 
has removed most of the high spatial-frequency noise, some 
traces of this noise are still visible in the image of the fil¬ 
tered data cube. This suggests that this filtering could be 
improved. One possible strategy is to use a Butterworth fil¬ 
ter, with n — 2, given by the following equation: 


H (u, v) = < 


1 + 


^) 2 +(^) 2 


( 2 ) 


which corresponds to the product between two identical el¬ 
liptical filters. The use of the Butterworth filter given by 
equation (2) is equivalent to apply two successive Butter¬ 
worth filterings with the same elliptical filter. If this process 
is performed with n = 2, no significant ring will be intro¬ 
duced in the filtered images, but the efficacy of the filtering 
will be improved. Again, we assume a = b in equation (2), 
which is equivalent to use a squared circular filter. It is im¬ 
portant to mention, however, that the cut-off frequency of 
a squared circular filter has to be a little higher than the 
cut-off frequency of a simple circular filter with the same 
value of n, in order to not compromise the PSF of the ob¬ 
servation. In the case of the data cube of the standard star 
observed with spaxel scale of 100 mas, our tests revealed 
that the most appropriate value for the cut-off frequency of 
the squared circular filter is 0.26 Ny. Fig. Illl shows images of 
an intermediate-wavelength interval of the Wo data cubes of 
this star before the Butterworth spatial filtering, after the 
Butterworth spatial filtering using a simple circular filter 
with n = 2 and cut-off frequency of 0.21 Ny, and after the 
Butterworth spatial filtering using a squared circular filter 
with n = 2 and cut-off frequency of 0.26 Ny. 

Fig. Illl shows that the Butterworth spatial filtering us¬ 
ing a squared circular filter with n — 2 removed a little more 
high spatial-frequency noise than the Butterworth spatial fil¬ 
tering using a simple circular filter with n = 2. Therefore, 
we conclude that the Butterworth spatial filtering using a 
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Figure 10. Moduli of the inverse Fourier transforms of the circular Butterworth filters, given by equation (1) (with a = b), with orders 
(a) n = 2, (b) n = 3, (c) n = 4, (d) n = 5 and (e) n = 6. All the filters were obtained using a cut-off frequency of 0.21 Ny. 
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Figure 11. Images of an intermediate-wavelength interval of the Wo data cube of HIP 048480 (with spaxel scale of 100 mas) (left) before 
the Butterworth spatial filtering, (centre) after the Butterworth spatial filtering using a simple circular filter with n = 2 and a cut-off 
frequency of 0.21 Ny and (right) after the Butterworth spatial filtering using a squared circular filter with n = 2 and cut-off frequency of 
0.26 Ny. Note that some traces of the high-frequency noise are still visible in the image from the data cube filtered with a simple circular 
filter, but are absent in the image from the data cube filtered with a squared circular filter. 


squared circular filter with n = 2 is the most appropriate 
to be applied to SINFONI data cubes. Although the previ¬ 
ous discussion was based on the results obtained with the 
data cube of the standard star observed with the fore-optics 
with spaxel scale of 100 mas, the conclusions are valid for 
observations taken with any fore-optics. 


The seeing conditions of the observing night and the 
effect of the AO (which does not have the same efficiency in 
all observations) may affect the value of the cut-off frequency 
to be used in the Butterworth spatial filtering. In the case of 
the fore-optics with an FOV of 8.0 arcsec (in which the AO 
is frequently not used), the values of the cut-off frequencies 
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Figure 12. Images of intermediate-wavelength intervals of the data cubes of (top) HIP 049220, (middle) HIP 048480 and (bottom) 
HIP 044245. The images at left represent the data cubes before the Butterworth spatial filtering. The images at centre represent the 
data cubes after the Butterworth spatial filtering with cut-off frequencies of (top) 0.29 Ny, (middle) 0.26 Ny and (bottom) 0.27 Ny. The 
images at right represent the data cubes corresponding to the difference between the non-filtered and filtered ones. 


will depend basically on the seeing conditions of the night 
(and also on the size of the spaxels of the instrument, which, 
as explained in Section 7, affect the PSF of the observation). 
In the case of the other two fore-optics, however, the cut-off 
frequencies will depend mainly on the efficiency of the AO 
(and, again, on the size of the spaxels of t he instrument). In 
our previous experiences jMenezesI 120121 1 . we have applied 
the Butterworth spatial filtering to SINFONI data cubes 
using cut-off frequencies between 0.20 and 0.35 Ny. 

Using a filter given by equation (2) and n — 2, we ap¬ 
plied the Butterworth spatial filtering to the original data 
cubes (without any wavelet decomposition) of the standard 
stars observed with spaxel scales of 25, 100 and 250 mas, 


with cut-off frequencies of 0.29, 0.26 and 0.27 Ny, respec¬ 
tively. Fig. 1121 shows images of intermediate-wavelength in¬ 
tervals of these data cubes, before and after the Butterworth 
spatial filtering, and also of the data cubes corresponding 
to the difference between the filtered and non-filtered data 
cubes of these objects. We can see that most of the high 
spatial-frequency noise (in the form of horizontal stripes) 
was removed from the images of the data cubes by the filter¬ 
ing procedure, without the introduction of detectable arte¬ 
facts. 

One important topic to be mentioned here is about the 
origin of the high spatial-frequency components in the im¬ 
ages of the data cubes with spaxel scales of 25, 100 and 250 
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Figure 13. Left: image of an intermediate-wavelength interval of the non-filtered data cube of HIP 048480 (with spaxel scale of 100 
mas), with spaxels of 50 mas. Centre: image of the same intermediate-wavelength interval of the filtered data cube of HIP 048480, with 
spaxels of 50 mas. The filtering was applied using a squared circular filter with n = 2 and a cut-off frequency of 0.55 Ny. Right: image 
of the same intermediate-wavelength interval of the data cube corresponding to the difference between the previous two. 


mas. As discussed in Paper I, such high-frequency compo¬ 
nents are not introduced by the interpolation we apply after 
the spatial re-sampling. We applied a Butterworth spatial fil¬ 
tering to the non re-sampled data cube of the star observed 
with spaxel scale of 100 mas (using a squared circular filter, 
with n=2), in order to determine if these high-frequency 
components were completely introduced by the spatial re¬ 
sampling or were, at least some of them, present in the non 
re-sampled data cubes. In this case, we verified that the most 
appropriate value for the cut-off frequency to be used in the 
Butterworth spatial filtering is 0.55 Ny. Fig.[l3]shows images 
of an intermediate-wavelength interval of the non re-sampled 
data cube of this star, before and after the Butterworth spa¬ 
tial filtering, and also of the data cube corresponding to the 
difference between the filtered and non-filtered ones. 

Fig. [13] reveals that the non re-sampled data cube of 
the star observed with spaxel scale of 100 mas has essen¬ 
tially the same high spatial-frequency components observed 
in the re-sampled data cube (Fig. 1121) . Therefore, we con¬ 
clude that, although the spatial re-sampling procedure has 
certainly introduced high spatial-frequency components in 
the images of this data cube (as the Fourier transforms in 
Fig- [7]show), most of them were already present in the orig¬ 
inal data cube. This discussion was based, only, on the data 
cube with spaxel scale of 100 mas; however, the results ob¬ 
tained with data cubes in the other fore-optics are entirely 
analogous. 


6 ‘INSTRUMENTAL FINGERPRINT’ 
REMOVAL 

After the Butterworth spatial filtering, the next step in our 
data treatment is the ‘instrumental fingerprint’ removal. Un¬ 
like the NIFS instrumental fingerprint we described in Pa¬ 
per I, which is relatively subtle and may not be relevant for 
most studies, the SINFONI instrumental fingerprint is sig¬ 
nificantly more intense and may affect many studies. This 
feature usually takes the form of a large horizontal stripe 
at the bottom of the images and shows a characteristic low- 
frequency spectral signature. One example of a study af¬ 
fected by the SINFONI instrumental fingerprint is the work 


made bv lNeumaver et all (120071 ) of data cubes of the nuclear 
region of the Centaurus A galaxy. These data were obtained 
in the H and K bands, using the fore-optics with an FOV 
of 3.2 arcsec. The [Fe II] images in that study were signifi¬ 
cantly affected by the fingerprint. The process for removing 
this feature is the same we use for removing instrumental 
fingerprints of NIFS data cubes (see Paper I ) and is b ased 
on the PC A Tomography technique dsteiner et al.ll2009l) . 

PCA Tomography involves the use of the Principal 
Component Analysis (PCA - Murtagh & Heck 1987; Fuku- 
naga 1990) technique on data cubes. In this application of 
PCA, the spectral pixels of the data cube are taken as the 
variables of a coordinate system and the spaxels are taken 
as the observables. PCA then creates a new system of un¬ 
correlated coordinates (called eigenvectors), which are linear 
combinations of the original coordinates (the wavelengths) 
and are orthogonal to each other. These eigenvectors are or¬ 
dered in a way that eigenvector El is associated with the 
highest fraction of the data variance, eigenvector E 2 is asso¬ 
ciated with the second highest fraction of the data variance 
and so on. Due to the way the eigenvectors are defined in the 
case of PCA Tomography, we call them eigenspectra. The 
projections of the observables (spaxels) on the eigenvectors 
are images, which we call tomograms. 

The process for removing the instrumental fingerprint 
of SINFONI data cubes involves the following steps: 

• removal of the spectral lines (both emission and absorp¬ 
tion lines) of the data cube, 

• application of PCA Tomography to the data cube with¬ 
out emission lines, 

• selection of the eigenvetors related to the instrumen¬ 
tal fingerprint (the eigenvectors related to the instrumental 
fingerprint must have the spectral signature of the finger¬ 
print and the corresponding tomograms must have the spa¬ 
tial morphology of a large horizontal stripe at the bottom 
of the images), 

• fit of splines to the eigenspectra related to the instru¬ 
mental fingerprint, 

• construction of a data cube using the splines fitted to 
the eigenspectra related to the instrumental fingerprint and 
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Figure 14. Tomogram and eigenspectrum corresponding to eigenvector E 4, obtained with PCA Tomography of the data cube of 
Centaurus A, with spaxel scale of 100 mas, after the removal of the spectral lines. The spline fitted to the eigenspectrum is shown in 
red. This eigenvector is dominated by the instrumental fingerprint, but the AGN emission can be detected in the tomogram. However, 
the effect of this emission is of second order. 
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Figure 15. Results of the subtraction of the image corresponding to the wavelength interval 2.1481 - 2.1530 /jm from the image 
corresponding to the wavelength interval 2.0988 - 2.1037 /im of the data cube of Centaurus A, with spaxel scale of 100 mas. The image 
at left was obtained before the instrumental fingerprint removal and the image at right was obtained after the fingerprint removal. 


the corresponding tomograms (the resulting data cube will 
contain only the fingerprint), 

• subtraction of the data cube obtained above from the 
original one. 

The discussion of all previous steps can be found in 
Paper I. In order to give an example of the efficacy of our 
method for instrumental fingerprint removal on SINFONI 
data cubes, we used a data cube of the nuclear region of 
Centaurus A, which was observed with SINFONI, in the 
K band, using the fore-optics with an F OV of 3.2 arcsec. 
This i s one of the data cubes analysed by iNeumaver et al.l 
(l2007h and it was treated with all the techniques described 
before. After the removal of the spectral lines, we applied 
the PCA Tomography on the resulting data cube. After a 
careful analysis, we verified that the obtained eigenvector 
E 4 reveals very clearly the instrumental fingerprint typically 
found in SINFONI data cubes. The eigenspectrum and the 
tomogram corresponding to eigenvector E 4 are shown in 
Fig. 1141 The spline fitted to the eigenspectrum is shown in 


red. Eigenvector E 4 explains ~ 7.0 x 10 -4 per cent of the 
variance of the data cube obtained after the removal of the 
spectral lines. 

Eigenspectrum E 4 in Fig. 1141 and the fitted spline, re¬ 
veal very clearly the spectral signature of the instrumental 
fingerprint. This spectral signature shows two main broad 
features: one around 2.15 /im and the other around 2.35 
/im. The feature at 2.15 /im is the most characteristic of the 
instrumental fingerprint of SINFONI data cubes obtained 
(with any fore-optics) in the K band. The feature at 2.35 
/im may also appear, although it is more diffuse in some 
cases. The small fraction of the data variance explained by 
eigenvector E 4 suggests that this instrumental fingerprint is 
not intense enough to compromise any analysis or measure¬ 
ment. However, as mentioned before, this ar tefact may affect 
consi derably the images of the data cubes (INeumaver et al.l 
l2007lf . In order to illustrate the effect of this fingerprint 
on a scientific analysis, we constructed images of the wave¬ 
length intervals 2.0988 - 2.1037 /rm and 2.1481 - 2.1530 /im 
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Figure 16. Tomograms and eigenspectra corresponding to the eigenvectors related to the instrumental fingerprints in the data cubes of 
Centaurus A, with spaxel scale of 250 mas, after the removal of the spectral lines, in the (top) J, (middle) H and (bottom) K bands. 
The splines fitted to the eigenspectra are shown in red. Although the eigenvectors are dominated by the instrumental fingerprint, the 
AGN emission can be seen in the tomogram obtained from the data cube in the K band and the shape of the diffuse stellar component 
is visible in all tomograms. However, as mentioned for Fig. 1141 the effects of this emission from the AGN and from the stellar component 
are of second order. 


of the data cube of Centaurus A. The second image was 
then subtracted from the first one. This procedure was ap¬ 
plied to the data cube before and after the instrumental 
fingerprint removal. The obtained images, shown in Fig. 1151 
reveal the considerable improvement provided by our cor¬ 
rection. The instrumental fingerprint can be easily seen in 
the non-corrected image, but essentially no fingerprint is de¬ 
tected in the corrected image. Since the fingerprint is associ¬ 
ated with the spectral continuum and not with the emission 


lines, analysis involving the continuum are usually more af¬ 
fected than other types of analysis. One important study 
of the spectral continuum that could be affected by the in¬ 
strumental fingerprint is the spectral synthesis. However, as 
mentioned before, in some cases, even studies involving the 
images of emission line s may be affected by the fingerprint 
dNeumaver et al.ll2007l ). 

In order to compare the observed instrumental finger¬ 
print in different spectral bands, we also used data cubes 
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Figure 17. Top: results of the subtraction of the image corresponding to the wavelength interval 1.2091 - 1.2156 /im from the image 
corresponding to the wavelength interval 1.1727 - 1.1792 /im (left) before and (right) after the instrumental fingerprint removal of the 
data cube of Centaurus A. Middle: results of the subtraction of the image corresponding to the wavelength interval 1.5905 - 1.6007 (im 
from the image corresponding to the wavelength interval 1.5480 - 1.5582 /im (left) before and (right) after the instrumental fingerprint 
removal of the data cube of Centaurus A. Bottom: results of the subtraction of the image corresponding to the wavelength interval 2.1368 
- 2.1454 fi m from the image corresponding to the wavelength interval 2.1098 - 2.1184 jin\ (left) before and (right) after the instrumental 
fingerprint removal of the data cube of Centaurus A. All data cubes were obtained with the fore-optics with spaxel scale of 250 mas. 


of the nuclear region of Centaurus A, observed, in the J, H 
and K bands, using the fore-optics with an FOV of 8.0 arc- 
sec. We chose this specific fore-optics because it was the only 
one with available observations in the three spectral bands of 
this object. We followed the procedure described above and 
applied PCA Tomography to these three data cubes, after 
the removal of the spectral lines. Fig. H6l shows the obtained 
eigenspectra and tomograms associated with the fingerprint 


in the data cubes. The splines fitted to the eigenspectra are 
shown in red. Eigenvectors E 3 (in the J band), E2, (in the 
H band) and E4 (in the K band) explain ~ 0.142, ~ 0.075 
and ~ 0.003 percent, respectively, of the data variances of 
the corresponding data cubes, after the removal of the spec¬ 
tral lines. We can see that the spatial morphology and the 
spectral signature of the fingerprint in the three spectral 
bands are actually very similar. The fingerprint has the ap- 
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pearance of a large horizontal stripe at the bottom of the 
images and the two broad features of the spectral signature 
are located around 1.15 and 1.22 /im, in the J band, around 
1.52 and 1.60 /im, in the H band, and, as mentioned before, 
around 2.15 and 2.30 /im, in the K band. As Figs. 14 and 
16 reveal, the characteristics of the instrumental fingerprint 
do not change with the fore-optics used in the observation. 
Again, in order to show the effect of the fingerprint on a sci¬ 
entific analysis of these data cubes, we constructed images 
of the wavelength intervals 1.1727 - 1.1792 /im and 1.2091 

- 1.2156 /im , in the J band, 1.5480 - 1.5582 /im and 1.5905 

- 1.6007 /im, in the H band, and 2.1098 - 2.1184 /im and 
2.1368 - 2.1454 /im, in the K band. Then, for each band, 
we subtracted the second wavelength interval from the first 
one. We applied this procedure before and after the instru¬ 
mental fingerprint removal. Fig. [lT] shows the results. It is 
easy to see the substantial improvements provided by our 
treatment procedure. 

The spectral characteristics of the eigenspectra in Figs. 
14 and 16, together with the spatial morphologies of the cor¬ 
responding tomograms, clearly indicate that these eigenvec¬ 
tors are mainly associated with the instr umental fingerp rint. 
However, in our previous experiences dMenezej [20121 1. we 
verified that eigenvectors dominated by a given phenomenon 
may also show contaminations by other phenomena. The to¬ 
mograms obtained from the data cubes of NGC 5643 in the 
K band, shown in Figs 14 and 16, for example, have struc¬ 
tures related to the active galactic nucleus (AGN) emission. 
The shape of the diffuse stellar component may also be seen 
in the tomograms from the data cubes in the J, H and K 
bands, shown in Fig. 1161 However, the effects of this emission 
from the AGN and from the stellar component are of sec¬ 
ond order, as the subtraction of the data cubes constructed 
using the eigenvectors associated with the fingerprint from 
the original data cubes removed the fingerprint but did not 
result in measurable changes in the observed emission. In 
other words, we can say that the eigenvectors in Figs 14 
and 16 are dominated by the instrumental fingerprint but 
also show contaminations by the central emission from the 
galaxy. Nevertheless, these contaminations did not interfere 
significantly in the process of instrumental fingerprint re¬ 
moval. 

The fractions of the data variances explained by the 
eigenvectors shown in Fig. [16] indicate that the influence of 
the fingerprint in the data cubes of Centaurus A decreases 
from the J band to the K band. In order to explain that, 
first of all, it is important to mention that, in our previous 
experiences {m eneze S 2012:), we have observed that objects 
with a more intense continuum emission are usually less af¬ 
fected by the instrumental fingerprint than objects with a 
weaker continuum. Therefore, since the continuum emission 
from the AGN in Centaurus A increases from the J band 
to the K band, the fact that the relevance of the finger¬ 
print decreases from the J band to the K band is actually 
expected. On the other hand, since the fingerprint is associ¬ 
ated with the spectral continuum and not with the emission 
lines, objects without continuum emission are not affected 
by this problem. In data cubes of point-like sources (like HIP 
049220, HIP 048480 and HIP 044245), most of the FOV does 
not show significant continuum emission. As a consequence, 
in such cases, the instrumental fingerprint normally appears 
in a subtler way or may not appear at all. This is the reason 


why we used data cubes of a galaxy to describe the process 
for instrumental fingerprint removal, which may not be nec¬ 
essary in cases in which the fingerprint is too subtle to be 
clearly detected. 

So far, we were unable to determine the origin of the 
instrumental fingerprint in SINFONI data cubes. Its charac¬ 
teristics do not seem to change with time. One possibility is 
that it is caused by some illu mination pattern in the instru¬ 
ment (iNeumaver et al.]|2007l l. However, another possibility 
is that the fingerprint has been introduced during the data 
reduc tion. Nevertheless, our previous experiences llMenezesI 
l2012h proved that our method is capable of removing this 
artefact in essentially all cases. 


7 RICHARDSON-LUCY DECONVOLUTION 

The last step in our sequence for the treatment of SIN- 
FONI data cubes is the Ri chardson-Lucy deconvolution 
(iRichardso nl l 1972|; |Lucy||l974| l. A deconvolution is an iter¬ 
ative process that reverts the effects of a convolution. The 
image of a celestial object observed from the Earth’s sur¬ 
face corresponds to the original image convolved with the 
PSF due to the Earth’s atmosphere and to the instrument. 
Therefore, in the astronomical case, a deconvolution is ap¬ 
plied in order to obtain the original image of an object, given 
the observed image and the PSF of the observation. 

We verified that, among all the deconvolution tech¬ 
niques that can be applied to astronomical images (Starck 
& Murtagh 2006), the Richardson-Lucy is the one that pro¬ 
vides the best results for SINFONI data cubes, as it does not 
introduce detectable artefacts in the images. For more de¬ 
tails about the Ric hardson-Lucy deconvolu tion, see Paper I. 
The use of the AO jDavies &: Kaspeiil2012l l usually provides 
a PSF containing a central diffraction spike, with a shape of 
an Airy function, and a surrounding halo, with a shape of a 
Lorentzian function. However, in the case of SINFONI data 
cubes, the PSFs obtained with the fore-optics with FOVs of 
0.8 and 3.2 arcsec can be described by the sum of a Gaus¬ 
sian and a Lorentzian function (equation 13 in Paper I). On 
the other hand, the PSFs of SINFONI data cubes obtained 
with the fore-optics with an FOV of 8.0 arcsec (in which the 
AO is usually not applied) are well described by a simple 
Gaussian function. 

Similarly to what we discussed in Paper I for NIFS data 
cubes, it is probable that the existence of a central Gaussian 
component in the PSFs of SINFONI data cubes (obtained 
with the two fore-optics with the highest spatial resolutions) 
is due to the original size of SINFONI spaxels. For an 8 
m telescope, the values of the full width at half-maximum 
(FWHM) of the Airy function at 1.25, 1.65 and 2.2 /im (the 
mean wavelengths of the J, H and K bands, respectively) 
are 0.032, 0.043 and 0.057 arcsec, respectively. For observa¬ 
tions taken with the fore-optics with spaxels of 0.125 arcsec 
x 0.25 arcsec, all these FWHM values are obviously smaller 
than one spaxel. So, it would not be possible to resolve the 
Airy profile, even if the AO was used (which is usually not 
the case for this fore-optics). One interesting point is that 
the spatial sampling of this fore-optics is not adequate to 
precisely resolve even a brightness profile with an FWHM 
of 0.5 arcsec (which is a typical value of the seeing in this 
spectral region). Therefore, although the final PSF of an ob- 
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Figure 18. Vertical brightness profiles of images of intermediate-wavelength intervals of the data cubes of (top) HIP 049220, (middle) HIP 
048480 and (bottom) HIP 044245. The red curves represent simple Gaussian fits and the green curves represent Gaussian + Lorentzian 
fits. 


servation taken with this fore-optics has a Gaussian shape, 
this is probably not a precise reproduction of the brightness 
profile of the PSF generated by the atmosphere. For obser¬ 
vations taken with the fore-optics with spaxels of 0.05 arcsec 
x 0.1 arcsec, the FWHM values are smaller than one of the 
dimensions of the SINFONI spaxels (0.1 arcsec) and very 
close to the other (0.05 arcsec). As a consequence, one can 
say that the spatial sampling is still not capable of resolving 
the Airy profile and this is the reason why a Gaussian pro¬ 
file is usually a good description for the central spike of the 
PSFs obtained with this fore-optics. For observations taken 
with the fore-optics with spaxels of 0.0125 arcsec x 0.025 
arcsec, the situation is a little different. The FWHM of the 
Airy profile in the J band is smaller than three spaxels along 
the horizontal axis and smaller than two spaxels along the 
vertical axis. Similarly, the FWHM of the Airy profile in the 


H band is smaller than four spaxels along the horizontal axis 
and smaller than two spaxels along the vertical axis. Finally, 
the FWHM of the Airy profile in the K band is smaller than 
five spaxels along the horizontal axis and smaller than three 
spaxels along the vertical axis. So, this fore-optics provides 
a better spatial sampling for re solving the Air y profile. How¬ 
ever, our previous experiences dMenezesll2012l ) revealed that, 
even in this case, the spatial sampling is usually not suffi¬ 
cient to distinguish between a Gaussian profile and an Airy 
profile and, consequently, the central spike can still be ade¬ 
quately described by a Gaussian function. One should note 
that the AO corrections work less perfectly in the J and H 
bands. As a consequence, in these two bands, it is difficult 
to reach the diffraction limit and the spaxel scales of 25 and 
100 mas are usually adequate to sample the obtained PSFs. 
Fig. [TS] shows the vertical brightness profiles of images of 
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Figure 19. Images of intermediate-wavelength intervals of the data cubes of (top) HIP 049220, (middle) HIP 048480 and (bottom) HIP 
044245. The images at left are from the data cubes before the Richardson-Lucy deconvolution and the ones at right are from the data 
cubes after the Richardson-Lucy deconvolution. 


intermediate-wavelength intervals of the data cubes of the 
standard stars observed with spaxel scales of 25, 100 and 
250 mas. Gaussian and Gaussian+Lorentzian fits are also 
shown. 

Fig. 1181 confirms that a Gaussian+Lorentzian fit is the 
most adequate to reproduce the vertical brightness profiles 
of the stars with spaxel scales of 25 and 100 mas and that 
a Gaussian fit reproduces, with good precision, the vertical 
brightness profile of the star with spaxel scale of 250 mas. 

Unlike the Airy function, the SINFONI spatial sampling 
is not wavelength dependent. Therefore, since the PSFs of 
the observations taken with the fore-optics with FOVs of 0.8 
and 3.2 arcsec depend significantly on the spatial sampling 


(which is not adequate to resolve the Airy profile), these 
PSFs remain approximately constant along the spectral axis 
of the data cubes. In the case of observations taken with the 
fore-optics with an FOV of 8.0 arcsec, as discussed before, al¬ 
though the corresponding PSFs have a Gaussian shape, the 
spatial sampling is not adequate to precisely reproduce the 
brightness profile of the PSF generated by the atmosphere. 
Therefore, the PSFs of observations taken with this fore¬ 
optics are also significantly affected by the spatial sampling 
and, consequently, remain approximately constant along the 
spectral axis. Considering this stability of the PSFs, we can 
say that a Richardson-Lucy deconvolution can be applied to 
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SINFONI data cubes (obtained with any fore-optics) using 
a constant PSF. 

The PSFs of SINFONI data cubes obtained with the 
fore-optics with FOVs of 0.8 and 3.2 arcsec are consider¬ 
ably difficult to reproduce, as they show a Gaussian and 
a Lorentzian component with a relative weight that varies 
from one observation to the other. Due to this complex¬ 
ity and the fact that the Richardson-Lucy deconvolution is 
very sensitive to the knowledge of the PSF, the best re¬ 
sults of the use of this procedure in these two fore-optics are 
achieved using real PSFs. As explained in Paper I, images 
of an intermediate-wavelength interval of data cubes of iso¬ 
lated stars or images of the broad wings of permitted lines 
in type 1 AGNs can be used as PSFs. If a real PSF is not 
available, one can construct a synthetic PSF as a sum of 
a Gaussian and a Lorentzian function. In the case of SIN¬ 
FONI data cubes obtained with the fore-optics with FOV 
of 8.0 arcsec, as explained before, the PSFs have a simple 
Gaussian shape. So, in these cases, it is relatively easy to 
construct synthetic PSFs and the Richardson-Lucy decon¬ 
volution performed with such PSFs provides results essen¬ 
tially as precise as the ones obtained using real PSFs. For 
all fore-optics, there are cases in which it is difficult to ob¬ 
tain an estimate of the PSF from the science data cube. In 
cases like that, it is possible to estimate the PSF from the 
data cube of the standard star used in the data reduction. 
This approach, however, should be performed with caution 
as the PSF of the data cube of the standard star may not 
be the same as the PSF of the data cube of the observed 
object, due to differences in the effects of the AO applied 
to these data cubes and also due to differences in the seeing 
during these observations. Nevertheless, in many situations, 
the differences between the PSFs of the science data cube 
and of the standard star data cube may be subtle and a 
reliable Richardson-Lucy deconvolution can be applied us¬ 
ing the PSF of the standard star data cube. When it is not 
possible to obtain any estimate of the PSF, the Richardson- 
Lucy deconvolution should not be applied. As established 
in Paper I, the Richardson-Lucy gives optimized results if a 
number of iterations between 6 and 10 is used. 

We applied the Richardson-Lucy deconvolution to the 
data cubes of the standard stars observed with spaxel 
scales of 25, 100 and 250 mas, using real PSFs (cor¬ 
responding to images of the data cubes in intermediate 
wavelengths) and 10 iterations. Fig. [T^] shows images of 
intermediate-wavelength intervals of the deconvolved and 
non-deconvolved data cubes of these objects. The Strehl ra¬ 
tios improved from 0.39 to 1.25 in the data cube with spaxel 
scale of 25 mas, from 0.04 to 0.10 in the data cube with 
spaxel scale of 100 mas and from 0.008 to 0.022 in the data 
cube with spaxel scale 250 mas. Without the Richardson- 
Lucy deconvolution, we usually obtain Strehl ratios in the 
ranges of 0.3 - 0.7 and 0.01 - 0.2 for data cubes in the fore¬ 
optics with FOVs of 0.8 and 3.2 arcsec, respectively. How¬ 
ever, with the Richardson-Lucy deconvolution, it is possible 
to obtain Strehl ratios as high as 0.3 for the fore-optics with 
FOV of 3.2 arcsec and e ven higher than 2 for the fore-optics 
with FOV of 0.8 arcsec jMavll201ll 'l. In the case of the fore¬ 
optics with a FOV of 8.0 arcsec, it is difficult to establish 
a range of values for the Strehl ratio, because these val¬ 
ues depend fundamentally on the seeing conditions of the 
night. The results obtained above indicate that the use of 


the Richardson-Lucy deconvolution in SINFONI data cubes 
provides a considerable improvement of the spatial resolu¬ 
tion, in any fore-optics. 

The fact that we obtain Strehl ratios higher than 1 in 
deconvolved data cubes in the fore-optics with the highest 
spatial resolution may be somewhat surprising, as, in any 
observation, this value should always be smaller or equal 
to 1. Indeed, such high values of the Strehl ratio are only 
possible here because of the use of deconvolution. Since this 
procedure is not limited by the diffraction laws as in a real 
observation, there is, in principle, no upper limit for the 
Strehl ratio of a deconvolved image._ 

In our previous experiences dMeneze i5m3) , we have 
never observed Strehl ratios higher than 0.7 in the fore¬ 
optics with an FOV of 0.8 arcsec and higher than 0.2 in 
the fore-optics with an FOV of 3.2 arcsec. These observed 
limits are probably related to the size of SINFONI spaxels 
and not only to limitations in the AO. Following the same 
strategy used for NIFS data cubes in Paper I, we calculated 
the Strehl ratios of images containing Airy functions, in the 
J, H and K bands, sampled in spaxels of 0.0125 arcsec x 
0.025 arcsec and 0.05 arcsec x 0.1 arcsec (the sizes of the 
spaxels in the two fore-optics with the highest spatial res¬ 
olutions). For the fore-optics with spaxels of 0.0125 arcsec 
x 0.025 arcsec, the calculated Strehl ratios are ~ 0.6, ~ 0.7 
and ~ 0.8, for the J, H and K bands, respectively. On the 
other hand, for the fore-optics with spaxels of 0.05 arcsec x 
0.1 arcsec, the calculated Strehl ratios are ~ 0.05, ~ 0.09 
and ~ 0.20, for the J, H and K bands, respectively. Based 
on these results, we can say that the observed limits in the 
Strehl ratios are probably a consequence of the size of SIN¬ 
FONI spaxels. 

8 A SCIENTIFIC EXAMPLE: NGC 5643 

NGC 5643 is an almost face-on (i ~ 27° ±5°) SAB(rs)c 
galaxy. Its nucleus shows a narrow emission-li ne spectrum 
(FWHM ~ 300 km s -1 ) of high ionization (ISandagelll978l L 
Based on these characteristics, this object was classified as 
a low luminosity Seyfert 2 by Phillips, Charles & Baldwin 
(1983). 

Several studies have shown that NGC 5643 seems to be 
an example of an object containing an obscured active nu¬ 
cleus with anisotropic escape of radiation. T his scenario is 
compa tible with the unified model for AGNs. I Morris et al. I 
(1198511 analysed velocity maps of optical emission lines, ob¬ 
tained with the imaging Fabri-Perot spectrometer TAURUS, 
and also flux maps in 6 cm and 20 cm, obtained with the 
Very Large Array. The authors observed an extended op¬ 
tical line emission with an elongated morphology, approxi¬ 
mately along the east-west direction (which coincides with 
the direction of the bar). They also found evidences that the 
velocity field is similar to the one of models of gas flux in a 
barred potential. The radio data revealed a nucleus coincid¬ 
ing with the optical nucleus and with the rotation centre, 
with two lobes (one at each side) along the same direction 
of the optical emission. Based on these observations, the au¬ 
thors proposed a simple model, in which the gas flows along 
the bar, forms a disc orthogonal to the major axis of the 
bar and collimates the radio emission, obscuring, also, the 
ionizing continuum. 
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Figure 20. Br 7 images from the (a) original and (b) treated data cubes of NGC 5643, and the corresponding RGB composite images 
from the (c) original and (d) treated data cubes of NGC 5643. The colours blue, green and red in panels (c) and (d) correspond to the 
velocity ranges —305 km s — 1 ^ V r ^ —14 km s' 1 , 0 km s —1 7 Vy 7 28 km s —1 and 42 km s~ 1 7 V T 332 km s~ 1 , respectively. 


ISchmitt. Storchi-Bergmann fc Baldwin! (11994 ”) analysed 
images and long-slit spectra, obtained with the Cerro Tololo 
Inter-American Observatory , and verified that the distri¬ 
butions of the emission-line ratios support the hypothe- 
sis that the io n izing continuum is collimated in a bicone. 
ISimpson et al.1 (1997?) analysed high resolution images of 
NGC 5643 obtained with the Hubble Space Telescope ( HST) 
and observed that the [O III]/Ha map has a well defined 
structure in the form of a ‘V’, which probably corresponds 
to the projection of a tridimensional cone. 

The data cube of NGC 5643 we analyse here is described 
in further detail in Section 2 and it was treated with all the 
procedures described before. [Hicks et ahl J2013! '). using stan¬ 
dard methodologies, analysed this data cube, together with 
the data cubes of other nine galaxies (all observed with SIN¬ 
FONI, in the K band, with the fore-optics with an FOV of 
8.0 arcsec), forming a sample containing five Seyfert galax¬ 
ies and five quiescent galaxies. The authors verified that the 
Seyferts in the sample have a more centrally concentrated 
nuclear stellar surface brightness, a lower stellar velocity dis¬ 
persion within a radius of 200 pc, an elevated H 2 1-0 S(l) 
luminosity out to a radius of at least 250 pc and a more 
centrally concentrated H 2 surface brightness. The authors 


interpreted these characteristics of the Seyfert galaxies as 
an indication for the existence of a thermally cold (in com¬ 
parison to the bulge) nuclear structure, composed of a gas 
reservoir and a relatively young stellar population, in these 
objects. This indicates that the fuelling of a Seyfert galaxy 
is associated with the formation of the dynamically cold nu¬ 
clear component. In this work, we are reanalysing the data 
cube of NGC 5643 with our data treatment methodology, 
showing different scientific aspects. It is important to em¬ 
phasize that none of the results obtained with our analysis 
of the data cube of NGC 5643 is in conflict with what was 
obtained bv lHicks et al.1 (120131 '). as the work made by these 
authors was focused on a comparison between Seyfert and 
inactive galaxies, not on the details of individual galaxies. 
Our analysis here is only about NGC 5643 and focused on 
the emission and kinematics of the ionized and molecular 
gas. 

We compared images from the treated and non-treated 
data cubes of NGC 5643. The non-treated data cube corre¬ 
sponds to the median of the five data cubes obtained after 
the data reduction. Only the correction of the DAR was ap¬ 
plied to each one of these five data cubes. Fig. [20] shows 
images of the Br 7 emission line (with the adjacent spec- 
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Figure 21. (a) [O III] image of the nuclear region of NGC 5643, obtained with WFPC2 of the HST, with the same FOV of the data 
cube analysed here, (b) The same image shown in (a), with a different look-up table, (c) The same image shown in (a) and (b), after a 
convolution with an estimate of the PSF of the data cube of NGC 5643. (d) Bry image from the treated data cube of NGC 5643. The 
main ionized-gas clouds are identified with different notations in (a)/(b), (c) and (d). 


tral continuum subtracted) from the treated and non-treated 
data cubes of NGC 5643. The corresponding RGB compos¬ 
ite images (based on the radial velocity intervals) are also 
shown. All the images show the existence of an extended 
emission approximately along the east-west direction. How¬ 
ever, only the images from the treated data cube allow a 
clear visualization of individual ionized-gas clouds, which 
are almost undetectable in the images from the non-treated 
data cube. The extended emission along the east-west direc¬ 
tion represents the narrow-line region (NLR), in the form of 
a bicone, of this galax y. This result i s com patible with the 
observations made by Morris et al. I (Il985l l and Schmitt et 


al. (1994). 

Fig. l21l shows an [O III] image with two different look-up 
tables of the nuclear region of NGC 5643, obtained with the 
Wide-Field Planetary Camera 2 (WFPC2) of the HST (this 
image was previously analysed by Simpson et al. 1997). Since 
the spatial resolution of the HST is considerably higher than 
the spatial resolution of the data cube of NGC 5643, it is 
difficult to establish a comparison between the [O III] image 
and the Br 7 image from the treated data cube. Therefore, 
in order to make that comparison, we convolved the [O III] 


image with an estimate of the PSF of the data cube. The 
result of this convolution and the Bry image from the treated 
data cube of NGC 5643 are also shown in Fig. [21] The main 
ionized-gas clouds are identified in all the images in Fig. 1211 
(and are compared in Table [2j. 

The first thing we can notice in Fig. [5T]is that the [O 
III] i mage reveals the exi stence of an one-sided ionization 
cone (ISimpson et al.lll997l ). as the ionized-gas clouds extend 
mainly to the east, but not to the west of the nucleus. Only 
one very weak ionized-gas cloud (21/i) is located at the east 
of the nucleus. On the other hand, the Br 7 image shows a 
bicone. The difference between these two images is probably 
due to the extinction by dust, which affects more the [O III] 
image than the Br 7 image. Indeed, Schmitt et al. (1994) con¬ 
cluded that there is an increase of the interstellar extinction 
to the west of the nucleus and a decrease of the interstellar 
extinction to the east. This behaviour is certainly compat¬ 
ible with the existence of an apparent one-sided ionization 
cone at the east of the nucleus, in the [O III] image. 

By analysing the positions of the ionized-gas clouds in 
Fig. 1211 we can see a good correspondence between the [O 
III] and the Br 7 images. There are, however, a few discrep- 
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Table 2. Comparison between the central ionized-gas clouds of NGC 5643 identified using SINFONI data and HST data. 


Identification of the cloud 
in the Br 7 image 

Identification of the cloud 
in the [O III] image convolved 
with an estimate of the PSF 
of the data cube of NGC 5643 

Identification of the cloud 
in the original [O III] image 

Radial velocity obtained 
using the Br 7 emission 
line (km s —1 ) 

B 

b 

5+6+7+8-1-9+10-1-11 

187 ± 18 

C 

c 

20 

184 + 5 

D 

d 

12 

-132 + 6 

E 

e ? 

13 ? 

44 + 9 

F 

/ 

14+15 

43 + 2 

G 

g+h ? 

16+17+18+19 ? 

-43 ± 26 

H 

i ? 

21 ? 

-111 ± 32 
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Figure 22. Br 7 emission line of the spectrum of the ionized- 
gas cloud C identified in the treated data cube of NGC 5643. A 
Gaussian fit is shown in red. 


ancies that should be discussed. Clouds g and h in the [O III] 
image, for example, appear as cloud G only in the Bry im¬ 
age. Clouds e and E appear to be associated with each other 
but their positions are a little different. The same happens 
to clouds i and H. These discrepancies may be related to the 
extinction by dust. However, it is important to mention that 
the [O III] emitting areas do not necessarily correspond to 
the Br 7 emitting areas. Therefore, the discrepancies between 
the [O III] and the Bry images in Fig.[2T]may represent real 
differences between the [O III] and the Bry emitting areas. 
Table [2] shows a comparison between the identifications of 
the ionized-gas clouds in all the images in Fig. [21] 

We extracted spectra from circular regions, with a ra¬ 
dius of 0.2 arcsec, centred on all the clouds of the NLR visible 
in the Bry image from the treated data cube. We then shifted 
all the spectra to the rest frame, using the task 1 dopcor’ 
(from the ‘noao’ package) with z = 0.003999 (NASA Ex- 
tragalactic Database - NED), and applied a correction due 
to the heliocentric velocity, using the task ‘rvcorrecV (from 
the ‘astutil’ package). After that, we fitted Gaussian func¬ 
tions to the Br 7 emission lines in the obtained spectra, in 
order to determine the radial velocity values of the ionized- 


gas clouds. This entire analysis was performed in iraf en¬ 
vironment. Fig. [22] shows the Gaussian fit applied to the 
Bry emission line of the spectrum of the ionized-gas cloud 
C. The profile of the emission line is well described by the 
Gaussian function and, therefore, we conclude that the fit 
provides a reliable value for the corresponding radial veloc¬ 
ity. The radial velocity values of the ionized-gas clouds are 
shown in Table [2] A velocity map of the Bry emission line 
would give more details about the kinematics of the ionized 
gas in the nuclear region of NGC 5643. However, we verified 
that the S/N ratio of the Bry emission line in the regions 
where most of the ionized-gas clouds are located is consid¬ 
erably low. As a consequence, the obtained velocity values 
are not sufficiently precise to construct a velocity map. This 
is the reason why we determined the velocity values of the 
Bry line for each ionized-gas cloud using spectra extracted 
from circular regions, centred on each cloud. 

Tabic [2] shows that the moduli of the radial velocities 
of the three ionized-gas clouds with the highest projected 
distances from the central AGN (clouds E, F and G) are 
considerably lower than the moduli of the radial velocities 
of the other clouds. We cannot perform a detailed analysis of 
the ionized-gas kinematics due to our limited FOV. However, 
the contrast between the radial velocities of clouds E, F and 
G and the others suggests the existence of an acceleration of 
the ionized gas near the AGN, followed by a deceleration due 
to the interaction with the interstellar medium. A similar 
behaviour was also observed in the NIFS data cube of the 
nuclear region of the Seyfert 1 galaxy NGC 4151 (see Paper 

I)- 

In order to analyse the molecular emission of H 2 at the 
surroundings of the AGN in NGC 5643, we constructed im¬ 
ages of the H 2 A 212 I 8 emission line (with the adjacent spec¬ 
tral continuum subtracted) from the treated and non-treated 
data cubes of this galaxy. Fig. l23l shows the obtained images 
and also RGB composite images based on the radial velocity 
intervals. Similarly to what was observed in the images of 
the Bry emission line in Fig. 1201 all the H 2 A 212 I 8 images 
in Fig. [23] show a clear extended emission; however, only 
in the images from the treated data cube we can identify 
individual molecular-gas clouds. Besides that, we can see 
that the highest velocity values of the molecular gas can be 
found approximately along the north-south direction, which 
is perpendicular to the direction of the ionization cone (east- 
west). The images in Fig. [23] may give the impression that 
part of the extended emission has been taken away by the 
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• 2-1012 - 2-1012 

Arcsec Arcsec 


Figure 23. H 2 A 212 I 8 images from the (a) original and (b) treated data cubes of NGC 5643, and the corresponding RGB composite 
images from the (c) original and (d) treated data cubes of NGC 5643. The colours blue, green and red in panels (c) and (d) correspond 
to the velocity ranges —345 km s 1 7 V r —48 km s 1 , —34 km s’ 1 ’ V r 7 0 km s —1 and 8 km s — 1 7 V r ■ 305 km s’ 1 , respectively. 


data treatment; however, this is just a consequence of the 
improvement of the spatial resolution due to the use of the 
Richardson-Lucy deconvolution. 

A model capable of explaining the observed morphol¬ 
ogy of the molecular gas in Fig. [53] involves the existence 
of a molecular torus/disc structure along the north-south 
direction. One possibility is that the gas flows along the 
torus/disc towards the nucleus, feeding the AGN. This 
torus/disc collimates the AGN emission along the ioniza¬ 
tion cone. This scenario i s com patible with the models pro¬ 
posed by iMorris et al. I (|l985l h Schmitt et al. (1994) and 
ISimpson et alJ (1997) and also with the unified model for 
AGNs. Fig. 1241 shows an RG composite image with the Br7 
emission in red and the H2A212I8 emission in green. This 
image makes it easier to visualize the proposed scenario in¬ 
volving a molecular torus/disc feeding the AGN and colli¬ 
mating the AGN emission. One important point that should 
be discussed here is the fact that Fig.[24]shows the existence 
of some molecular gas in areas corresponding to the ioniza¬ 
tion cone of NGC 5643. Considering that most of the H2 
molecules should be destroyed along the ionization cone, it 
is probable that the superposition of emission from ionized 


and molecular gas we are seeing in Fig. [24] is caused by the 
projection on the plane of the sky and the corresponding 
emitting areas are actually not coplanar. A simple test for 
that hypothesis would be to compare the velocity values 
of the ionized-gas clouds and molecular-gas clouds that are 
apparently superposed. If the values are compatible, then a 
reasonable assumption is that the superposed areas are asso¬ 
ciated with each other. On the other hand, if the values are 
not compatible, then it is very unlikely that the superposed 
areas are related to each other and the apparent superposi¬ 
tions are probably just caused by the projection on the plane 
of the sky. In order to perform this test, we extracted the 
spectra from circular regions, with a radius of 0.2 arcsec, cen¬ 
tred on each molecular-gas cloud detected in the H2A212I8 
image from the treated data cube of NGC 5643. The identi¬ 
fications of the molecular-gas clouds in the H2A212I8 image 
are shown in Fig. 1251 After that, we repeated the procedure 
used with the spectra extracted from the ionized-gas clouds: 
we shifted all the spectra to the rest frame (using the task 
‘ dopcor ’ from the ‘noao’ package, with 2 = 0.003999) and 
applied a correction due to the heliocentric velocity (using 
the task 1 rvcorrect ’ from the ‘astutil’ package). Then, we 
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Table 3. Radial velocities calculated for the main molecular-gas clouds identified in the H2A21218 image from the treated data cube of 
NGC 5643. 


Identification of the cloud 
in the H2A212I8 
image 

Radial velocity obtained 
using the H2A212I8 
emission line (km s —1 ) 

B’ 

89 ±2 

C’ 

42 ± 2 

D’ 

-57 ±2 

E’ 

-47 ± 2 

F’ 

-17 ±5 

G’ 

28 ±3 

H’ 

-11 ±4 




Figure 24. RG composite image from the treated data cube 
of NGC 5643, with the Br7 emission in red and the H2A212I8 
emission in green. 

fitted Gaussian functions to the H 2 A 212 I 8 emission lines, 
in order to determine the radial velocities of the molecular- 
gas clouds. Fig. [7S] shows the Gaussian fit applied to the 
H 2 A 212 I 8 of the spectrum of the molecular-gas cloud C’. 
The radial velocity value provided by this fit is essentially 
as precise as the one provided by the Gaussian fit applied 
to the Br 7 emission line in Fig. 1221 The radial velocities 
obtained for the molecular-gas clouds are shown in Table [3] 
Table [3] confirms that the molecular-gas clouds with the 
highest velocity values are located approximately along the 
north-south direction. Figs. 21, 24 and 25 reveal that there 
is a partial superposition between clouds C and F’ and be¬ 
tween clouds F and H’. No more significant superpositions 
can be observed. However, Tables 2 and 3 show that the 
velocity values of clouds C and F are not compatible (at 
3a level) with the velocity values of clouds F’ and H\ re¬ 
spectively. This indicates that, as expected, the molecular 
emission of H 2 is not coplanar with the ionized-gas emission 
along the ionization cone of the NGC 5643. 

9 SUMMARY AND CONCLUSIONS 

We presented a treatment procedure for VLT/SINFONI 
data cubes, which includes the following steps: correction 


Figure 25. H2A212I8 image from the treated data cube of NGC 
5643, with the main molecular-gas clouds identified. 



Wavelength (tun) 

Figure 26. H2A212I8 emission line of the spectrum of the 
molecular-gas cloud C’ identified in the treated data cube of NGC 
5643. A Gaussian fit is shown in red. 
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of the DAR; spatial re-sampling, which provides a better 
visualization of the contours of the spatial structures; But- 
terworth spatial filtering, to remove high spatial-frequency 
components from the images of the data cubes; instrumen¬ 
tal fingerprint removal; Richardson-Lucy deconvolution, to 
improve the spatial resolution of the data cubes. 

The main conclusions related to the steps of the data 
treatment procedure are as follows. 

• The DAR effect may be significant for SINFONI data 
cubes obtained with the fore-optics with FOVs of 0.8 and 3.2 
arcsec, even if only one spectral band is considered. In the 
case of the fore-optics with FOV of 8.0 arcsec, the correc¬ 
tion of the DAR may not be necessary. Both the theoretical 
and the practical approaches are very effective in removing 
the DAR, although the practical approach is a little more 
precise. 

• The spatial re-sampling followed by an interpolation of 
the values provides a clearer visualization of the contours 
of the spatial structures. We usually apply the spatial re¬ 
sampling to SINFONI data cubes (observed with any fore¬ 
optics) in order to double the number of spaxels along the 
horizontal and vertical axis. 

• A Butterworth spatial filtering using a filter corre¬ 
sponding to the product between two identical circular fil¬ 
ters removes most of the high spatial-frequency components 
(mainly in the form of horizontal stripes) from SINFONI 
data cubes, without compromising the images. 

• An instrumental fingerprint in the form of a large hori¬ 
zontal stripe at the bottom of the images and with a charac¬ 
teristic spectral signature is usually found in SINFONI data 
cubes. Our method involving the PCA Tomography tech¬ 
nique is very effective in removing this artefact from the 
data cubes. 

• With the Richardson-Lucy deconvolution, it is possible 
to obtain Strehl ratios as high as 0.3 for the fore-optics with 
an FOV of 3.2 arcsec and even higher than 2 for the fore¬ 
optics with an FOV of 0.8 arcsec. The PSFs of SINFONI 
data cubes correspond to a simple Gaussian function for 
observations obtained with the fore-optics with an FOV of 
8.0 arcsec. On the other hand, the PSFs are composed of a 
Gaussian and a Lorentzian component for observations ob¬ 
tained with the fore-optics with FOVs of 0.8 and 3.2 arcsec. 

In order to give a scientific example of the efficacy of 
our data treatment procedure, we analysed a data cube of 
the nuclear region of the galaxy NGC 5643, obtained with 
the fore-optics with an FOV of 8.0 arcsec, and compared the 
results obtained with and without the treatment procedure. 
The main conclusions of this analysis are the following. 

• An extended emission can be easily seen in the images 
of the Br 7 emission line from the treated and non-treated 
data cubes of NGC 5643. However, only the image from the 
treated data cube allows a clear visualization of individual 
ionized-gas clouds, which are almost impossible to visualize 
in the image from the non-treated data cube. 

• Although we cannot perform a detailed analysis of the 
ionized-gas kinematics (due to our limited FOV), the radial 
velocities determined for the ionized-gas clouds suggest the 
existence of an acceleration of the ionized gas near the AGN, 
followed by a deceleration due to the interaction with the 
interstellar medium. 


• An extended molecular emission can also be seen in 
the images of the H 2 A 212 I 8 emission line from the treated 
and non-treated data cubes of NGC 5643. However, only 
the image from the treated data cube, again, allows a clear 
visualization of individual molecular-gas clouds. 

• The morphology of the H 2 molecular emission can be 
explained by the existence of a molecular torus/disc struc¬ 
ture along the north-south direction. One hypothesis is that 
the gas flows along the torus/disc towards the nucleus, feed¬ 
ing the AGN. This torus/disc collimates the AGN emission 
along the ionization co ne. This scenario is c ompatible with 
the models proposed bvlMo rris et af~l (1 1985ll , Schmitt et al. 
(1994) and ISimpson et ah (1997) and also with the unified 
model for AGNs. 


Considering the benefits provided by our data treat¬ 
ment procedure, we conclude that this methodology may 
improve the quality of the analysis of SINFONi data cubes. 
The scripts of all procedures included in our method can be 
found at http://www.astro.iag.usp.br/~PCAtomography. 
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